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Aggregated Demand Response Modelling for Future 

Grid Scenarios 
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Abstract 

With the increased penetration of intermittent renewable energy sources (RESs) in future grids (FGs), balancing between 
supply and demand will become more dependent on demand response (DR) and energy storage. Thus, FG feasibility studies 
will need to consider DR for modelling nett future demand. Against this backdrop, this paper proposes a demand model which 
integrates the aggregated effect of DR in a simplified representation of the effect of market/dispatch processes aiming at minimising 
the overall cost of supplying electrical energy. The conventional demand model in the optimisation formulation is augmented by 
including the aggregated effect of numerous users equipped with rooftop photovoltaic (PV)-storage systems. The proposed model 
is suited for system studies at higher voltage levels in which users are assumed to be price anticipators. As a case study, the 
effect of the demand model is studied on the load profile, balancing and loadability of the Australian National Electricity Market 
in 2020 with the increased penetration of RESs. The results are compared with the demand model in which users are assumed to 
be price takers. 

Index Terms 

Aggregated demand modelling, battery storage, demand response, future grids, photovoltaic generation, renewable energy 
sources. 

Nomenclature 

n £ N Supplier n and set of N suppliers. 
to € M Load aggregator m and set of M load aggregators. 
h £ H Time slot h and set of II slots in the horizon. 
ri £ R Node i in the system and set of R Nodes. 

i,j Indices. 

a, /? Incremental changes for and calculation. 

Ah Time step interval. 

Pq (h) Generated active power by supplier n in slot h. 

Cq(Pq) Cost function of supplier n. 
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Bid of supplier n for generating Pq . 

Maximum power limit of supplier n. 

Minimum power limit of supplier n. 

Total energy requirement of aggregator m over a horizon. 

Aggregated flexible power demand of aggregator m in slot h. 

Maximum flexible demand limit of aggregator m. 

Minimum flexible demand limit of aggregator m. 

Aggregated inflexible demand of aggregator to in slot h. 

Aggregated demand of price-responsive users communicating with aggregator to before reshaping in slot h. 
Aggregated PV generation of price-responsive users communicating with aggregator m in slot h. 

Aggregated nett demand of aggregator to in slot h. 

Battery storage power for aggregator to in slot h. 

Maximum battery charging rate for aggregator m. This is a limiting variable to ensure that the total storage capacity 
is not exceeded. 

Battery storage state of charge (SOC) of price-responsive users communicating with aggregator to,. 

Maximum battery storage SOC limit of price-responsive users communicating with aggregator to. 

Minimum battery storage SOC limit of price-responsive users communicating with aggregator to. 

Total battery energy loss of price-responsive users communicating with aggregator m over a horizon. 

Transferred power by line from r, to r ? in slot h. 

Maximum line power limit from r, to r y . 

Minimum line power limit from r t to r :/ . 

Susceptance of line between r, and /y. 

Voltage angle at node r, in slot h. 


I. Introduction 


F UTURE grid (FG) feasibility studies have demonstrated that relying on high penetration of diverse renewable energy 
sources (RESs) is possible assuming enough backup generation and/or utility storage are available to keep the network 
in balance Q-@. A preliminary study by the University of Melbourne Energy Research Institute has proposed a zero-carbon 
electrical grid for Australia in 2020 0- The University of New South Wales researchers have analysed the viability of 100 % 
RES scenarios considering a copper plate model for the Australian National Electricity Market (NEM) 0. 0. They have 
suggested 100 % RESs electricity in the NEM, at the current reliability standard, would be technologically feasible. Also, the 
least-cost mix of 100 % RESs scenario has been determined for the future of the NEM. Similarly, the least-cost mix of high 
penetration of diverse RESs and conventional generation has been determined for the future of the PJM, California and New 
Zealand networks in 0-@, respectively. 


However, these studies have only focused on simple balancing by using a simplified grid model such as the copper plate 
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model. On the other hand, the penetration of distributed generation (DG) has been increasing significantly in recent years, 
and greater penetration of battery storage is anticipated in power systems 0-G3- In particular, global installed capacity of 
rooftop photovoltaic (PV) has increased from approximately 4 GW in 2003 to nearly 128 GW in 2013 mainly due to electricity 
price increases, government incentives and also worldwide drop of PV capital cost 0 , 0 - In Australia, installed capacity of 
rooftop PV (which is mostly installed by residential and commercial customers) has grown from approximately 0.8 GW in 
2011 to over 4 GW in 2014 [ 10 ]. Recent studies have demonstrated that users equipped with PV-battery systems will reach 
retail price parity in the foreseeable future in the USA grids and the NEM 0-0- In light of these developments, a question 
arises how to model the aggregated nett demand (including DG, storage and demand response (DR)) to study FG scenarios. 

While the effect of DR is neglected in most of the existing FG feasibility studies [Tj-0, it is considered in few studies 
mainly through two different ways: 

Implicit modelling : DR is considered implicitly, but it is not reflected into the loads. For instance in 0, the effect of DR is 
considered through improving the capacity credit value for intermittent RESs (i.e. intermittency of RESs is decreased). However, 
due to the significant effect of loads on performance and stability of power systems, it can be expected that incorporating DR 
explicitly into the load models will affect the results of FG feasibility studies significantly. 

Explicit modelling : In few recent studies (HMH)’ the aggregated effect of DR is reflected into the conventional demand 
models. Including the effect of DR into the demand models requires allowing for the interaction between demand and supply 
sides in some ways. This is mainly done through three different approaches. Firstly, in some studies, the supply-side is modelled 
physically while price-responsive users are not represented physically 0, frg. in 0, the effect of flexible loads is analysed 
on reserve markets. That study presumes flexible load to be represented precisely by a tank model. Also, the reserve market 


is too simplified and physical constraints of the electrical grids (e.g. line limits) are not considered. In 1151, flexible demand 
is represented via a price-elasticity matrix. The elasticities are a measure of the change in demand in response to a change in 
the electricity price, and are typically obtained from the analysis of historical data. Secondly, there are some other studies that 
model demand-side technologies physically while the supply-side is represented through the electricity price profile 0 - That 
model assumes user to be price takers, i.e. the effect of user actions is not considered in the electricity price. Thirdly, in few 
recent studies, both demand-side technologies and supply-side are modelled physically 03, 03- This approach necessitates 
the need for integrated simulations in which both supply and demand sides need to be optimised jointly. As discussed in 03’ 
this approach can provide more realistic results. In fl8) , the aggregated charging management approaches for plug-in electrical 
vehicles (PEVs) is integrated into the market clearing process. The market process, however, is too simplified and physical 
constraints of the electrical grids are not considered. 

Although those models have proven their merits, a generic modelling approach is still required which can be used for any 
granularity level in the grid (e.g. from a city, a state or even the whole network). Furthermore, for studying FG scenarios of 
uptake of the various demand-side technologies, it is necessary to have a model which can easily integrate them. Against this 
backdrop, in this paper, we make a further step by proposing a model that integrates the aggregated effect of DR in a simplified 
representation of the effect of market/dispatch processes. The proposed optimisation formulation aims at minimising the overall 
cost of supplying electrical energy in which the conventional demand model is augmented by including the aggregated effect of 
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numerous price anticipating users equipped with PV-storage systems. Due to the price anticipating assumption, load aggregators 
are considered implicitly. As a case study, the effect of the demand model on the load profile, balancing and loadability of the 
NEM is studied using a 14-generator model ED- The results are also compared with the demand model in m in which users 
are assumed to be price takers. Eight scenarios are analysed with one business as usual (BAU) and seven different DR scenarios 
with renewable integration. For the BAU Scenario in 2020, the electricity supply is dominated by coal, gas, hydro, and biomass; 
and in the Renewable Scenarios, some of the conventional coal generators in Queensland and South Australia are replaced 
with wind farms (WFs) and concentrated solar plants (CSPs) with storage, as suggested in 0> ED to meet Australia’s RES 
target ]20) . Simulation results show that increasing the penetration of DR with price anticipating users, improves loadability 
with the increased intermittent supply penetration in the grid. 

The remainder of the paper is organised as follows: Section II proposes the aggregated demand model considering DR. Section 
III describes the test-bed assumptions and modelling. Section IV describes simulation scenarios, and discusses simulation results. 
Finally, Section V gives conclusions. 

II. Aggregated Demand Model Considering DR 

Aggregated load models are commonly used in system studies to reflect the combined effect of numerous physical loads ED’ 
|22 ;|. These can be inspired by physical devices, e.g. using a large induction motor to represent all the motors connected, or by 
data-driven approaches. Conventional aggregate load models only account for the accumulated effect of numerous independent 
load changes and some relatively minor control actions. Including the effect of DR requires allowing for the interaction between 
demand and supply sides in some way, e.g. price signals. In this study, we propose a model that integrates the aggregated effect 
of DR in a simplified representation of the effect of market/dispatch processes inspired by the traditional unit commitment 
problem. The proposed optimisation problem aims at minimising the overall cost of supplying electrical energy in which the 
conventional demand model is augmented by including the aggregated effect of numerous users equipped with newer demand 
technologies. The demand model consists of two parts: (i) a fixed electricity demand profile (inflexible demand), and (ii) flexible 
demand equipped with demand-side technologies (we consider a large homogeneous population of residential and commercial 
PV-storage systems, but the model allows for an easy integration of other demand technologies as well). Also, the demand 
model considers the following assumptions: 

Assumption 1 : Users are assumed to be price anticipators, i.e. the effect of users actions is considered on the electricity price 
by the load aggregators. Due to the price anticipating assumption, load aggregators are considered implicitly. We assume that 
users are incentivised by the load aggregators through a proper price signal derived using mechanism design [ |23| . 

Assumption 2: Aggregators do not aim to change the total energy consumption of the users, but instead to systematically 
manage and shift it. Furthermore, all load aggregators are connected to not only the power grid but also to a communication 
infrastructure to enable two-way communication with the users. 

Assumption 3: Users have smart meters equipped with smart home energy management systems (SHEMSs). SHEMS 
implements an algorithm to schedule energy resources and storage, and so achieves DR. 
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A. Optimisation model 

The optimisation model aims at minimising the overall electricity cost. The objective function of the model can be written 
as: 

H N 

min EE^( p G(M), a) 

h =1 n—1 

where, each decision horizon for the model (i.e. 24-hour period) is divided into one hour time-steps, giving a total of H = 24 
time-steps; denote a particular time-step by h , subject to the following constraints: 

1) Power generation limit: Generation of each supplier is a decision variable, and is constrained between the minimum and 
the maximum power limits as follows: 

P G <P G (h)<P G v/i, n, (2) 

2 ) Flexible demand, storage and PV: The following set of equations augment the conventional demand model by including 
the aggregated effect of numerous users equipped with PV-storage systems. 

p mm,m < < pmax.m ^ ^ ( 3 a ) 

H 

YsIPlfW - P£(h)]Ah = E m + Vm, (3b) 

h=l 

Flexible demand of each load aggregator, /j!", is a decision variable which reflects the aggregated effect of DR, and is 
constrained between the minimum and the maximum limits in (3a). The overall energy balance over time horizon H is given 
by (3b). Fig. □ shows a simple illustration of the demand profile. As it can be seen in Fig. [TJi, the aggregated nett demand of 
each load aggregator, P j’J, is equal to the sum of inflexible and flexible demands, i.e. P£p(/i) = P™(h) + P™(h) . When p™ n m 
is equal or more than zero ( P"' < Pj^), the model represent a situation where price-responsive users would not send power 
back to the grid. This situation might happen in the future if feed-in tariffs (FiTs) are much less than the retail tariffs paid by the 
users. On the other hand, if there will be a reasonable incentive for users in the future to send power back to the grid, P™ m,m 
can be relaxed to a value below zero (P^g < !'{") to reflect such a situation. The flexibility of loads is due to battery storage 
which is modelled implicitly by considering the upper limit of the flexible load power as P™ ax,m (/i) = -fjj'cha" 1 Note 

that P^ 1 ^™ is a limiting variable that ensures that the total storage capacity is not exceeded, and it does not represent a physical 
property of a particular battery technology. This variable is determined heuristically, as explained later in this section. The battery 
storage power, Pg 1 , is a byproduct of the optimisation problem which can be calculated as P^(h) = P™(h) — P™(h) — Pw{h). 
Fig. & shows the aggregated flexible demand equipped with PV-storage systems. As it is demonstrated in Fig. |TJ?, P" : '{h) is 
the aggregated demand of price-responsive users before utilising PV-storage systems. 

The aggregated flexible demand is determined by each aggregator in (3b) in a way that the total energy requirement for that 
aggregator remains constant ( Assumption 2) considering aggregated PV generation and also battery storage energy loss. Total 
energy requirement of aggregator m over the horizon can be written as E m = Y^h^iiPiTW + P™(h))Ah. E^" 1 in (3b) 
guarantees that battery round-trip efficiency is taken into account. As it is shown in Fig. [TJ), part of the required energy for 
price responsive users is provided by PV. So, the rest of the energy has to be supplied from the grid (i.e. + A”' — /ij"j. 
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Fig. 1. Simple illustration of the demand profile (a) aggregated flexible demand, inflexible demand and new nett demand profile, and (b) aggregated price- 
responsive demand and PV generation. 


Price-responsive users use their battery storage to shift their consumption to cheaper time slots to minimise the overall cost 
of supplying electrical energy (i.e. E™ + E™ — E™ can be spread out over the horizon due to enough storage capacity of 
price-responsive users in a way that the overall electricity cost minimises, the red area in Fig. [TJi). In other words, the total 
flexible demand energy over the horizon is equal to the energy which has to be supplied from the grid plus battery storage 
energy loss, i.e. £f=i P™(h)Ah = £f = 1 (P u m (/t) - Pft(h))Ah + = E? + E J* - E? + E»’™. 

3) Demand supply balance: Ignoring the losses in the grid, the power balancing equation can be written as (4a): 


E p s( h ) - E = E p l" Ti w v/i >^ e 


(h) = B ritr .(d ri (h) - 5 rj (h)) 


P, 


■ ^ / (ti) < P i ’ j 

L,min — 1/ — r L,max 


< p 


V/i, (ri,rj) G R , 
V/i, (ri,rj) G R, 


(4a) 

(4b) 

(4c) 


where, the power transferred between different nodes in the system is given by (4b), and is constrained by the line limits in 
(4c). 

To solve the above optimisation problem, aggregated demand, PV power, maximum battery storage charging rate and 
corresponding battery storage energy loss are required for each load aggregator. The next subsection describes how the maximum 
battery storage charging rate and corresponding energy loss can be calculated for the model. 


B. Modelling battery storage 

In the proposed optimisation formulation, battery storage is modelled implicitly and its state of charge (SOC) is not a 
decision variable. However, it is important to consider battery SOC limits to make sure that the available storage capacity 
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in the grid is not exceeded. is a limit parameter which ensures that the total storage capacity is not exceeded. It is 

determined using a heuristic search in Algorithm [T] 


Algorithm 1 : The proposed algorithm to determine P^^ 71 

1: Set 7 = 1, 

2: Initialize P~ = 0, 

3: Run Algorithm 2, 

4: while is within limits (7) do 
5: Solve the model (1-4). 

6: Calculate BQ q C (6). 

7: Z ^ — Z T 1 , 

„ n max.m , j-jmax.m , 

8: P B,cha,i P B,cha,(z-l) + 

9: Run Algorithm 2, 

10: end while 

11: Step 3: Return P" (l _ iy 


The battery energy loss corresponding to -P^ha? ' s also verihed based on a heuristic search in Algorithm |2j Based on the 
Algorithm [2j the optimisation formulation (1-4) is solved using the aggregated demand and PV power, suppliers’ data, battery 
storage charging rate (from Algorithm [lj and initial value for battery energy loss. Then, Ef 0 ”' can be calculated based on the 
storage actions as follows: 

H 

^i B or = (l-^)E P B m W ifP£'(h)> o, (5) 


h=l 


If the difference between the battery energy loss calculated from the equation (5) and the battery energy loss from the j th 
iteration, is less than the error, e, search stops for that horizon. Otherwise, the battery energy loss value in the j th 

iteration changes in a small step (/3) and the optimisation formulation solves again. The search continues until the difference 
between those two is less than the error. 


Algorithm 2 : The proposed algorithm to verify E^™ corresponding to l\ 


B,cha,i 


Set j = 1, 


Initialize E®’™. = 0, 
while > a do 

Solve the model (1-4). 
Calculate (5). 

j j + 1, 

pB,m pB ,m , n 

^lossj -^loss.O'-l) ” l ” 

end while 

Step 4: Return E*£.. 


After battery energy loss calculation from Algorithm [2] the optimisation formulation can be solved in Algorithm |T] for 
determining P %*To verify battery storage SOC, B£q C , which is a byproduct of the optimisation problem needs 

to be calculated: 


B 


m 

SOC 


(1 )=B 


t,m 

SOC’ 


(6a) 


B?oc(h) = B^ocih - 1) + P^{h - 1) Vft > 1, 


(6b) 










and compared with the battery storage SOC limits: 


nmin,m 

n soc 


< B™ c (h) < B™? 


Mh > 1, 


(7) 


While B^ c is within the limits, increases in a small step (a), corresponding to -P™haT calculates from 

Algorithm [2] and the optimisation formulation solves again. This procedure repeats until B$q C violates the limit. From a step 
before violation, -P™^" can be obtained. 


III. The Australian NEM Model 

A 14-generator model of the NEM, which was originally proposed for small signal stability studies jT9j, is used as the 
test-bed. The schematic diagram of the 14-generator model of the NEM is shown in Fig. [2] Areas 1 to 5 represent Snowy Hydro 
(SH), New South Wales (NSW), Victoria (VIC), Queensland (QLD) and South Australia (SA), respectively. The Australian 
Electricity Market Operator (AEMO) has proposed 16 zones for the NEM to capture differences in generation technology 
capabilities, costs, weather and so on in the future 0 - In order to extract data for the proposed model and generators in 
2020, the 14-generator model is matched with the 16 zones, as shown in Fig. [2] 

The generator technologies and the test-bed assumptions in this study follow Reference 0- In power system studies, 
generators’ cost functions are often modelled as quadratic. In this study, it is assumed that conventional generators submit their 
block bidding curves, while RESs are assumed to bid at zero cost. So, the cost function of the conventional generators can be 
approximated by a piece-wise linear function as: 

Ca{PS(h)) = Pa(PS{h))PS(h) Vn, (8) 

where, pg(^ c ( h )) are suppliers’ bids. Generators’ bids mimic the bidding behaviour of the generators in the NEM and use 
the predicted fuel price, thermal efficiency, and variable O&M costs 2020 0- 0- Table |T] lists bids of the suppliers (pi, p‘2 
and p3) in the 14-generator model in 2020. 

Each region of the 14 generator model (i.e. QLD, NSW, VIC and SA) is considered as a node for the demand model. 
Consequently, interstate lines are considered explicitly in the demand model. However, balancing and loadability studies are 
done using the 14 generator model shown in Fig. [2] The dispatch results from the market are used for balancing and loadability 
studies in DIgSILENT Power Factory. In the balancing studies, if supply cannot meet the demand, the hour is recorded as 
unserved hour. However, if available generation exceeds demand (i.e. due to high generation of RESs), the surplus power is 
recorded as dumped energy and that hour is marked as a dumped hour. 

IV. Simulation Scenarios and Results 

The effect of price anticipating and price taking assumptions, and different DR penetrations on the load profile, balancing 
and loadability of the NEM in 2020 with the increased penetration of RESs in the grid is demonstrated in this section. 
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Fig. 2. 14-generator model of the NEM. 


TABLE I 

THE GENERATORS’ BIDS IN 2020 


Gen. 

Type 

AEMO 

zone 

Pi 

($/MWh) 

P 2 

($/MWh) 

p3 

($/MWh) 

BPS2 

Coal 

NNS 

28.45 

42.66 

56.90 

EPS_2 

GT 

CAN 

69.20 

346.0 

692.0 

MPS_2 

Coal 

SWNSW 

27.43 

41.15 

54.86 

VPS_2 

Coal 

NCEN 

26.40 

39.60 

52.80 

LPS_3 

Biomass 

MEL 

39.50 

59.25 

79.00 

YPS_3 

Coal 

LV 

21.88 

32.82 

43.76 

CPS_4 

Coal 

CQ 

26.14 

39.21 

52.28 

GPS_4 

Coal 

CQ 

26.14 

39.21 

52.28 

SPS_4 

Coal 

NQ 

32.74 

49.11 

65.48 

TPS_4 

GT 

SWQ 

73.84 

369.2 

738.4 

NPS_5 

Coal 

NSA 

30.89 

46.34 

61.78 

PPS_5 

Coal 

SESA 

30.89 

46.34 

61.78 

TPS_5 

GT 

ADE 

69.20 

346.0 

692.0 
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TABLE II 

THE AGGREGATED STORAGE AND PV CAPACITIES FOR EACH REGION OF THE NEM FOR DIFFERENT UPTAKE SCENARIOS 


Region 

Scenario 

j^max,m 

&SOC ~^SOC 

(GWh) 

-pmax ,m 
B,cha 

(GW) 

PV capacity 
(GW) 


Low 

0.4-4.3 

0.44 

1.3 

QLD 

Medium 

0.6-6.4 

0.63 

1.9 


High 

0.9-8.5 

0.83 

2.6 


Low 

0.7-6.7 

0.62 

2.0 

NSW 

Medium 

1.0-10.1 

1.01 

3.0 


High 

1.4-13.5 

1.34 

4.1 


Low 

0.5-5.0 

0.47 

1.5 

VIC 

Medium 

0.8-7.5 

0.72 

2.3 


High 

1.0-10.0 

0.95 

3.0 


Low 

0.1-1.2 

0.10 

0.3 

SA 

Medium 

0.2-1.7 

0.16 

0.5 


High 

0.2-2.3 

0.22 

0.7 


A. Simulation scenarios 

Eight scenarios are analysed with one BAU and seven different DR scenarios with renewable integration. The Renewable 
Scenarios are analysed with the conventional load, and different levels of DR in the proposed demand model (Section II) and 
the demand model in GD- For the BAU Scenario, combinations of coal, gas, hydro, and biomass are considered for the NEM 
to supply the load in 2020. Then, some of the conventional coal generators in QLD and SA are replaced with CSPs together 
with storage and WFs, respectively to meet Australia’s RES target. Displacement of conventional generators in the Renewable 
Scenarios and chosen capacities for RESs are inspired by studies in q, d). NPS_5 in SA is replaced with a WF with the 
capacity of 3 GW using NSA data d). Also, SPS_4 and GPS_4 in QLD are replaced with two CSPs with the capacity of 
4.5 GW each and using NQ and CQ data dl’ respectively. It was found that delaying CSP output by 12 hours minimises the 
unserved and dumped energy. The RESs serve about 20 % of the total demand in the Renewable Scenarios. 

In this study, hourly demand and PV power for the proposed demand model are obtained from the AEMO predications 
for 2020 dJ- Also, DR is considered for the residential and commercial customers, i.e. 60% of the total system load in the 
NEM in 2020 |TT) . The industrial customers are left unaffected. Furthermore, the percentage of the residential and commercial 
customers with PV are considered 20 %, 30 % and 40 % for low, medium and high uptake scenarios, respectively. Table [11] 
shows the aggregated storage and PV capacities for each region of the NEM for different uptake scenarios. The chosen PV 
capacities for different uptake scenarios are inspired by the AEMO study d). Moreover, the chosen storage capacities roughly 
correspond to typical PV and storage capacity for a household in Australia. Using the algorithms in Section II, is also 

calculated for different DR scenarios, and is reported in Table [II] 


B. Accuracy of the load model 

Fig. I demonstrates the results of the demand model in NSW for low DR Scenario during the 30 th and 31 st of January 
2020. As it can be seen in Fig. [ 3 J 1 , flexible demand, Pp SW , increases during cheap hours in comparison with p^sw because 
price-responsive users use their battery storage to store cheaper electricity (Fig. [ 3 J 5 ). However, flexible demand decreases during 
peak hours due to PV generation, Tpy w , and also battery storage discharge (Fig. Hr). Also, Fig. 0 clearly illustrates that 
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Fig. 3. (a) Aggregated demand of price-responsive users before reshaping, P^ sw , aggregated flexible demand, P^ sw , and aggregated PV generation, P^ w , 
(b) battery charging, discharging and SOC in NSW during 30 th and 31 st of January 2020. 


Eti P™(h)&h = EtiiPvW Pft(h))Ah + E^. Furthermore, Fig. §, shows that with appropriate selection of the 
battery SOC remains within the bounds. This shows the efficacy of the proposed model.lt should be noted that all the 
results in this section assume that price anticipating users do not sent power back to the grid {P™ < Iff). 


C. Comparison of the load profile with different approaches 

Figs. & and compare the effect of price anticipating and price taking assumptions, and also different levels of DR on 
the load profile resulting from solving the proposed model and the demand model in |[l6| in NSW during the 25 th to 28 th of 
May 2020, respectively. As it is shown in the Figs. [4ji and [4j>, in both approaches, the users shift their consumption (using 
PV-storage system) from expensive hours to cheaper ones to utilise zero cost electricity produced by RESs. In the proposed 
model, the effect of the users’ action on the electricity price signal is considered by the load aggregators. As a result, the load 
profile with price anticipating assumption (Fig. |4|t) is smoother than the load profile with price taking assumption (Fig. [djt). In 
the latter approach the effect of users’ action is not considered on the electricity price signal. So, price-responsive users may 
shift their consumption to cheaper time slots all-together, which results in secondary load peaks due to load synchronisation, 
as it can be seen in Fig. |4j5. 

In the next subsection, the effect of the both demand models on balancing and loadability of the NEM in 2020 are studied 
with the increased penetration of RESs. 


D. Balancing and loadability results 


1) BAU Scenario: In the BAU Scenario, a big portion of the demand is supplied by coal-fired power plants and the peak 
loads are met with backup supply (i.e. gas turbines (GTs)). The results for spilled energy and hours, supplied electrical energy 


from GTs and average loadability for all the scenarios over the simulated year are summarised in Table III For the loadability 


calculation, we checked different load/generation increase scenarios. Among them, we showed the results for one of the critical 
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Fig. 4. Demand profile in NSW during 25 th to 28 th of May 2020 with (a) price anticipating and (b) price taking assumptions. 


ones in Table III in which all loads in QLD are assumed to increase uniformly in small steps with constant power factor. Also, 
it is assumed that all the generators in QLD are scheduled with the same participation factor to pick up the system loads. The 
loadability is computed for each hour until a step before power flow divergence. 


2) Renewable Scenarios: In the rest of the paper. Renewable Scenarios with the conventional load, the load model in 1161 
and the proposed model are called CL, PTDR and PADR, respectively. Also, PTDR and PADR scenarios with low, medium and 


high DR penetrations are called PTDR1-PTDR3 and PADR1-PADR3 respectively in Table III Unserved hours for all scenarios 
are zero. Comparing the BAU Scenario and the Renewable Scenario with the conventional load (CL), it can be seen that with 
the increased penetration of RESs, the loadability is reduced from 27.13 GW to 22.17 GW.Also, the required electrical energy 
from the backup generation (i.e. GTs) is increased from 18.73 TWto 18.77 TW. 

Compared to the CL Scenario (i.e. zero DR penetration), higher penetrations of DR (PADR and PTDR) improve the 
performance and loadability, and reduce the required energy from backup supply. High DR Scenario with price anticipating 
assumption (PADR3) has the highest loadability followed by the medium and low DR Scenarios. The loadability is increased 
from 22.17 GW for the CL Scenario to 26.24 GW for the PADR3 Scenario, which implies a considerable improvement in 
loadability. However, with price taking assumption, medium DR Scenario (PTDR2) has the highest loadability (i.e. 25.53 GW). 
Increasing the penetration of DR with price taking assumption beyond a certain point (from PTDR2 to PTDR3), fails to 
improve loadability further. This is because of secondary peaks for high DR scenario which deteriorates loadability compared 
to lower DR penetrations. Furthermore, in both approaches, high DR Scenario has the lowest spilled energy and hours followed 
by medium and low DR Scenarios. Also, the required energy from the backup generation for DR scenarios with price taking 
assumption is slightly higher than DR scenarios with price anticipating assumption. This happens as the price taker users may 
shift their consumption to the cheaper time slots all-together, and therefore some backup supply might be required to keep the 
network in balance. 


Fig. [5] shows the the effect of different DR scenarios on the load profile when users are price anticipators in one of the 
critical summer (i.e. 07 th to 10 th of February) and winter peaks (i.e. 19 th to 22 th of July) for the NEM. The balancing results 
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TABLE III 

BALANCING AND LOADABILITY RESULTS FOR ALL THE SCENARIOS 


Scenarios 

Spilled energy 
(TWh) 

Spilled hours 

(%) 

GT energy 
(TWh) 

Loadability 

(GW) 

BAU 

- 

- 

18.73 

27.13 

CL 

0.71 

13.65 

18.77 

22.17 

PADR1 

0.66 

13.03 

18.28 

24.82 

PADR2 

0.61 

12.67 

17.78 

25.86 

PADR3 

0.54 

11.71 

17.04 

26.24 

PTDR1 

0.66 

13.03 

18.33 

23.78 

PTDR2 

0.61 

12.67 

17.85 

25.53 

PTDR3 

0.54 

11.71 

17.28 

24.41 


for the high DR scenario during those peak hours are demonstrated in Fig. [6] As shown in Fig. [6] in the critical summer 
days the wind is not strong enough and the output of WF is low, and in critical winter days the solar exposure reduces and 
CSP output decreases as well. However, due to enough response capacity from DR (Fig. [5J and backup supply, balancing and 
loadability under these worst-case conditions are maintained. Figs. [5] and [6] show that DR can help balance fluctuating RESs 
and demand in real-time even under the worst-case conditions. 


V. Conclusion 

This paper proposes a demand model which integrates the aggregated effect of DR in a simplified representation of the 
effect of market/dispatch processes. The model is intended to be used for system studies at transmission levels in which users 
are assumed to be price anticipators. In the proposed optimisation formulation, the conventional demand model is augmented 
by including the aggregated effect of numerous price anticipating users equipped with PV-storage systems. As a case study, 
the effect of the demand model on performance and loadability of the NEM in 2020 with the increased penetration of RESs 
is studied using a 14-generator model. Also, the results are compared with the demand model in which users are assumed to 
be price takers. 

Simulation results show that with the increased penetration of RESs and no DR loadability is reduced. With DR, however, 
loadability is improved and the required backup supply is reduced. Increasing the penetration of DR with price anticipating 
assumption improves loadability. However, when users are price takers, increasing the penetration of DR beyond a certain 
point does not necessarily improve the loadability and might even deteriorate it. This is due to load synchronisation of price 
taking users, which might results in secondary peaks. 

The future research aims at expanding the idea in this paper by using bi-level optimisation to avoid implicit modelling 
of battery storage and heuristic search. Also, the effect of the aggregated demand model including DR will be studied on 
performance and stability of FG scenarios. 
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